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Abstract 

We present the result of a calculation for the first even moment of the non-singlet four- 
loop anomalous dimension of Wilson twist-2 operators in QCD with full color and flavor 
structures. 



Calculation of anomalous dimensions of the Wilson twist- 2 operators is one of the part 
of operator product expansion for the structure functions in the framework of perturbative 
Quantum Chromodynamics (QCD). At the present time such calculations are performed 
up to three-loop order [H El [3J HI [5J El [7] , while other part of operator product expansion, 
the coefficient functions, which are known in the same order (HI [91 [TOj, [11], demand the 
four-loop anomalous dimensions. 

In this paper we present the result of calculations for the first even moment of the non- 
singlet anomalous dimension of Wilson twist-2 operators at fourth order in perturbative 
QCD. Similar result can be found in Ref. [12], but our result contains full color and flavor 
structures and the calculations are performed with a different method^. 

The moments of the structure functions are expressed through the parameters of 
the following operator product expansion of the T-product of electromagnetic currents: 



r„„ = i I d l z e""T{J„(z)J,(0)} 
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+ singlet contributions + higher twists , (1) 

where it is usual to use the following notation 

_ J?_ _ a, ,„x 
3 ~ 16tt2 _ 4vr { ' 

for the QCD strong coupling constant. The sum in eq. ([TJ runs over the standard set of 
the spin-iV, twist-2 irreducible (i.e. symmetrical and traceless in indices . . /ijv) flavor 
non-singlet quark operators: 

0«.{w.»mat} = ^A a 7 {m ^ M2 • • • V^ty, a = 1,2,..., 8, (3) 

where D ft are the covariant derivatives, A a are the generators of the flavor group SU(rif) 
and C% N {Q 2 1 [i 2 , a s ) are the corresponding coefficient functions. 

The non-singlet moments of the structure functions are expressed through operator 
product expansion ([T}) in the following form: 



M k , N = I d*x N - 2 Fr en (x,Q 2 ) 

Ck,N ( ~Tl a s J [^AT,proton(/ i2 ) — ^V,ncutron(/ i2 )] (4) 

n \ i / 



1 Note, that there is all-loop prediction for the 0(1/Nf) contribution to the non-singlet anomalous 
dimension of twist-2 operators in QCD |13j . 
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where A^nucieon is the spin-averaged nucleoli matrix elements of the operator: 

(p, nucleon|O a '^ 1 - w >|nucleon,p) = p^ 1 . . .p MJv} A^ nuclcon (/i 2 ) . (5) 

Application of the renormalization group technique gives for the coefficient functions 
the following standard expression: 



C a k , N a s (^ = CS>(1, a s (Q 2 )) exp (- £® } da' s ^ J . (6) 



)2 \ / ra 3 (Q 2 ) 

'a s (M 2 ) 

The anomalous dimensions 77V in eq. (jSJ) are defined as 



o d \ogZ N _ V> (n) ( n+ l) 



" n=0 

and renormalized operators and bare ones are connected as follows: 

(O a '^-^) R = (Zn)- 1 [p a ^-^) B . (8) 

The four-loop approximation for the /3-function in QCD in the MS-scheme was obtained 
in Refs. [HICE]: 

P(a s ) = -/3 a 2 s - ia 3 a - p 2 at - (3 3 a 5 s + O (af) , (9) 
/3 = ll-^n / , (10) 

ft = 102 - jn f , (11) 
n 2857 5033 325 , 

/3 2 = n f + n 2 f , (12) 

2 18 54 V ; 

„ / 149753 \ / 1078361 6508 \ 

A = (_ + 3564 Cs j - (^g- + Cs j », 

/50065 6472 \ 1093 

V 162 81 S 7 f 729 f K J 

The perturbative expansion in a s for the coefficient functions is 

C k , N (l, a.) = flfj, + fl&a. + + 4> s 3 + O (aj) , (14) 

where the CallanGross relation gives -B^ ^ = for all N and the standard deep inelastic 

normalization [8] of the coefficient functions implies B% N = 1. 

The next-next-next-to-leading (NNNL) approximations for the non-singlet moments 
after renormalization group improvement are 

M 2 , N (Q 2 ) = a^' /A) (b<% + B { 2 %a s + B^ N a 2 s + B^ N a^E(a s )A N (fi 2 ) (15) 



with 



E(a, 



o 



a s at 
1 + l3l El + W 



El 



where Ei = 7^ /?o — and f° r l ne NNNL approximation one should keep the four 

leading orders in the product of the power series for the coefficient functions with the 
series E(a s ). So, for the NNNL approximation to the non-singlet moments M 2j at we need 
to known the 3-loop coefficients B\ N for the coefficient functions 6*2,^(1, a s ) and the 4-loop 

coefficients 7^ for the anomalous dimensions 7jv(a s ). 

For the first even moment (N = 2) non-singlet operator ([3]) has the following form 



S. - E 2 M + E x ( {131 - /3 /5 2 ) /3 2 + E 2 (3l) - E 2 ^(3 3 + E 3 ft 
,(°) 



, (16) 
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where D = 4 — 2e is space-time dimension. 

The calculation of the anomalous dimension of such operators can be performed in a 
usual way through the computation of the Green's function with the operator insertion, 
which have the following general form in momentum space (see [T5]): 



,(2) 
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where p is the momentum flowing through the external quark legs. To determine different 
components we use the following projectors (see [16J) 
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Really, to find the anomalous dimension of the operator Cns we should compute only 
E^nsO )- -^ total number of four-loop diagrams is 12816. As in our previous work [T71 
ESI Ml US EI] all calculations were performed with FORM [22], using FORM package 
COLOR [23J for evaluation of the color traces. For the dealing with a huge number of 
diagrams we use a program DIANA [21] , which call QGRAF [25] to generate all diagrams. 
For evaluation of Feynman integrals we used the method from Refs. [26, 27] and our 
own implementation of the Laporta's algorithm [2H] in the form of the MATHEMATICA 
package BAMBA with the master integrals from Ref. 
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For the renormalization we need the three-loop renormalization constant for the oper- 
ator insertion g^jj\ a 'j^A v ^ with two quarks and one gluon legs, which can be obtained 
order by order in a usual way from the renormalization of operator as 



7 71/271/27 



(21) 



where Za, Z g and are the renormalization constants for gluon filed A^, coupling con- 
stant and quark correspondingly. 
Our final result is 
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(22) 
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where (see Ref. [15] ) 



d abcd d abcd N 4 _ gjyrt + lg 



N A ~ 96iV 2 

d abcd d abcd N C (N 2 C + 6) 
iVA ~ 8 ' 

d abcd d abcd jV 2 (iV 2 + 36) 

iV^ ~ 24 
and for the color group SU(N C ) the more simple Casimir operators are: 



(23) 
(24) 
(25) 



1 /V 2 — 1 

?> = -, CV = ^-, C A = N C , N A = N 2 c -l. (26) 

The obtained result coincides with the existing result for the first even moment of the four- 
loop non-singlet anomalous dimension from Ref. [12] if we substitute the explicit expression 
for all Casimir operators for QCD with three active quarks (i.e. for the gauge group SU(3) 
with rif = 3). Moreover, the part of our result, which is proportional to (n/) l_1 a*, coincide 
with the prediction from Ref. [13J, while the non-planar part was calculated by us in 
Ref. [21]. 

The substitution of the colour factors with N c = 3 into eq. (T2"2"|) gives the following result 
for QCD case 
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In conclusion we give the explicit results for the different number of active quarks: 

7ns'°° p (2, n/ = 3) = 3.55556 a s + 38.84774 a s 2 + 448.07162 a s 3 + 6532.13656 a s 4 , (28) 

7ns' oop (2, n f = 4) = 3.55556 a, + 35.68724 a s 2 + 306.02989 a, 3 + 3679.66906 a s 4 , (29) 

7 Ns k)op (2, n/ = 5) = 3.55556 a s + 32.52675 a s 2 + 161.53001 a s 3 + 1108.56696 a, 4 , (30) 

7^ oop (2, n f = 6) = 3.55556 a s + 29.36626 a s 2 + 14.571953 a s 3 - 1169.71912 a/ . (31) 
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It is interesting to compare our result with the predictions [30p, coming from the Pade 
resummation, which for 7^ 3 g(2,n^ = 4) gives 2629 or 2557 depending on the resummations 
procedure. Note, that in four-loop order new colour structures (|23p -( j25|) appear, which can 
disimprove resummation. So, we give below our result for 7ns (2) with the contributions 
from different colour structures: 

7 <g (2, n f = 3) = 4626.76262 + 1932.76417 - 27.39022 d F J , (32) 

7^(2,n/ = 4) = 1783.42519 + 1932.76417 - 36.52029 d$f , (33) 

7^(2,71/ = 5) = -778.54684+ 1932.76417 </- 45.65037 </, (34) 

7 <g(2, n f = 6) = -3047.70285 + 1932.764174/ ~ 54.78044 df/ , (35) 

where dff and are the contributions coming from (|23|) and f[2~4|) correspondingly. 
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